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Abstract 

We present a detailed physical analysis of the dynamical glass- 
jamming transition which occurs for the so called Knight models re- 
cently introduced and analyzed in a joint work with D.S.Fisher [15]. 
Furthermore, we review some of our previous works on Kinetically 
Constrained Models. 

The Knights models correspond to a new class of kinetically con- 
strained models which provide the first example of finite dimensional 
models with an ideal glass-jamming transition. This is due to the un- 
derlying percolation transition of particles which are mutually blocked 
by the constraints. This jamming percolation has unconventional fea- 
tures: it is discontinuous (i.e. the percolating cluster is compact at 
the transition) and the typical size of the clusters diverges faster than 
any power law when p y pc- These properties give rise for Knight 
models to an ergodicity breaking transition at pc'- at and above pc a 
finite fraction of the system is frozen. In turn, this finite jump in the 
density of frozen sites leads to a two step relaxation for dynamic cor- 
relations in the unjammed phase, analogous to that of glass forming 
liquids. Also, due to the faster than power law divergence of the dy- 
namical correlation length, relaxation times diverge in a way similar 
to the Vogel-Fulcher law. 
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1 Introduction 



The formation of amorphous sohds as glasses and granular media, i.e the 
glass and jamming transitions, are still unsettled and fascinating questions 
in condensed matter physics despite all the works that have been devoted to 
this subject. These phenomena occur in a variety of systems which, even if 
microscopically very different, share common features. Among the others we 
recall supercooled liquids, colloidal suspensions and non-thermal jamming 
systems (e.g. vibrated granular materials) [1-3]. Basic glassy properties 
include a dramatic slowing down of dynamics when a proper external pa- 
rameter is tuned (e.g. temperature is lowered for liquids) and the occurrence 
of a complicated relaxation: non exponential and spatially heterogeneous [4]. 
When relaxation times become longer than experimental scales, equilibrium 
can no more be achieved and the systems freeze into an amorphous phase. 

Even the basic issues in understanding these phenomena remain unsolved. 
In particular, it is not settled whether the dynamical arrest is due to the prox- 
imity of a phase transition and whether this putative phase transition is a 
static or purely dynamical one. However experiments make it clear that, if 
an ideal glass transition occurs, it should have an unconventional behavior 
with mixed first and second order features. On the one hand, the divergence 
of relaxation times and the fact that both entropy and internal energy seem 
continuous (for molecular liquids) is indicative of a second order transition. 
On the other hand, there is a discontinuous order parameter: the infinite 
time limit of the Fourier transform of the density-density correlation (anal- 
ogous to Edwards- Anderson parameter for spin glasses) has a finite jump at 
the transition. This corresponds to the fact that the modulation of the mi- 
croscopic density profile of the glass does not appear continuously from the 
fiat liquid profile (we generally refer for all systems to glass and liquid phase 
meaning the regime before and after the freezing into an amorphous solid). 
Besides these mixed first/second order properties, another unconventional 
feature, compared to usual phase transitions, concerns the scaling of relax- 
ation times. Different functional forms have been proposed in the literature. 
A very popular and successful representation is the Vogel-Fulcher law. This 
suggests that the logarithm of the relaxation time diverges as the inverse of 
the distance from the transition, i.e. logr ~ 1/T — Tq, for molecular liquids. 
Finally, one of the most puzzling features is the absence of any experimental 
evidence of a static diverging correlation length. In particular, the dramatic 
slowing down of dynamics does not seem to be due to an increasing long 
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range order as for e.g. ferromagnetic transition: typical glass configurations 
are not very different from instantaneous configurations of the liquid. An 
enormous amount of theoretical approaches have been proposed in the last 
fifty years to describe these physical phenomena. However, if one takes the a 
priori assumption that a real glass or jamming phase transition takes place 
at a finite temperature (at a density less than the close-packed one for the 
jamming transition) then the number of scenarii reduces drastically. Leav- 
ing aside very phenomenological ones, there remain only two: either a very 
subtle static transition occurs (random first order scenario [5]) or the transi- 
tion is purely dynamical. The first possibility invokes [5] a thermodynamic 
transition similar to the one occurring for the so called discontinuous spin 
glass models, e.g. REM and p-spin models [6]. This is a one step replica 
symmetry breaking transition with a discontinuous (Edwards- Anderson) or- 
der parameter but no discontinuity for energy and entropy. This approach 
has been successful in explaining and predicting many physical phenomena 
related to the glass transition [7]. Although mean field models upon which it 
is based are very well understood by now, some of the finite dimensional pre- 
dictions are still semi-phenomenological and need further theoretical works 
to be put on a firm and solid basis, see [8-11] for recent works in this di- 
rection. The possibility of a purely dynamical glass transition have been 
mainly investigated through the so called Kinetically Constrained Models 
(KCM) ( [12] — [41] and references therein). These are stochastic lattice 
gases based on the ansatz that glass or jamming transitions are due to effec- 
tive geometrical constraints on the rearrangements of the atoms or molecules 
generated close to the transition. Static correlations beyond those present 
in dense liquids are assumed to play no role. Recently, in a joint work with 
D.S.Fisher [15], we have shown that indeed some KCMs display a purely 
dynamical transition on finite dimensional lattices with the basic features of 
the previously described glass-jamming transition. 

Note that these two scenarii are not necessarily in contrast because it 
might well be that, despite similar behaviors, the glass-jamming transition 
in molecular liquids, colloids and granular media are of different origin. The 
thermodynamic scenario might apply to the glass transition of molecular 
liquids whereas jamming transitions might correspond to purely dynamic 
phase transitions. 

Here we give an extended explanation of the result in [15] and we review 
some of our previous works (and some others when needed) on KCM. We 
focus in particular on the different behaviors and tools needed to deal with 
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different choices of the stochastic dynamics. This paper is not intended to 
be an extensive review on KCMs: for all the results previous to our works we 
refer to the reviews [12, 13] and references therein, while some more recent 
results for kinetically constrained spin models can be found in [20-23]. 

The outline of the paper is the following. In Section |21 we introduce the 
models distinguishing those with Glauber (KCSM) and Kawasaki (KCLG) 
dynamics and further dividing both classes into cooperative and non cooper- 
ative models. In Section El we review the different types of dynamical arrest 
that occur and motivate the analysis of KCMs from a theoretical point of 
view. In Section |3] we present tools and results for non cooperative mod- 
els. In Section ini we describe the jamming transition that cooperative KCMs 
display on Bethe lattices. In Section [HI we discuss the class of cooperative 
models (Knight models) introduced in [15,37] proving that they display an 
ideal glass transition. In particular, we show that this transition is related 
to a new type of percolation transition, which we call jamming percolation. 
In this Section we also explain which are the tools one needs in general to 
study cooperative KCMs, using Knight models as an example. Finally, we 
present our conclusions in Section [7| 

2 Kinetically Constrained Models (KCM) 

KCMs are stochastic lattice gases with hard core exclusion, that is on each 
site there can be one or zero particle. In other words, a configuration on a 
lattice A is given by the set of two-valued occupation variables on each site 
X G A: r]x = {1,0}. These represent occupied and empty sites (particles and 
vacancies), respectively. The dynamics is given by a continuous time Markov 
process which consists of a sequence of jumps for models with conserva- 
tive (Kawasaki) dynamics and birth/death for models with non conservative 
(Glauber) dynamics. The former are also known as Kinetically Constrained 
Lattice Gases (KCLG), the latter as Kinetically Constrained Spin Models 
(KCSM) or facilitated spin models (77^ = {1,0} are interpreted as spin up 
and spin down). For all the models we consider, dynamics satisfies detailed 
balance w.r.t. Bernoulli product measure, fip , at density p. Thus, there are 
no static interactions beyond hard core and an equilibrium transition cannot 
occur. However, in order for a move to be allowed, it is not enough to verify 
the hard core constraint. Indeed the jump or birth/death rates are non zero 
only if the configuration satisfies some additional local constraints, hence the 
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name kinetically constrained. These mimic the geometric constraints on the 
possible rearrangements in physical systems, which could be at the root of 
the dynamical arrest [14,17]. For example, in a highly dense liquid, numer- 
ical simulations [38] have shown that a molecule is caged by its neighbors 
and it cannot move on a substantial distance unless the neighbors move and 
the "cage is opened" . Experimental evidence of such cage effect for colloidal 
suspensions and granular materials have indeed been detected [3,39]. 

Numerical simulations show that, for proper choices of the constraints, 
KCM display glassy features including stretched exponential relaxation, super- 
Arrhenius slowing down and dynamical heterogeneities [13]. Therefore, de- 
spite their simplified and discrete character they could capture the key in- 
gredients of glassy dynamics, at least on proper time scales. Indeed, several 
works have been recently devoted to understanding the mechanism which 
induces these glassy properties and evaluating the typical time/length scales 
involved. 

2.1 Kinetically Constrained Spin Models (KCSM) 

KCSMs are endowed with a non conservative dynamics: each site changes 
from occupied to empty and from empty to occupied with rate (1 — p)fx{v) 
and pfx{v)y respectively. The value of fxijf)^ which encodes the kinetic con- 
straints, does not depend on rj^. Thus detailed balance is satisfied with 
Bernoulli product measure, /ip, at density p. 

KCSMs can be divided into two classes : non-cooperative and cooperative 
ones. For the former it is (for the latter it is not) possible to construct an 
allowed path which completely empties any configuration provided a proper 
finite cluster of vacancies is present somewhere. We call this cluster a defect. 
Among non cooperative models we recall the Fredrickson-Andersen [14] one 
spin facilitated (FAl/). For FAl/ a move in x is allowed only if at least 
one of the nearest neighbors is empty: fx{v) = 1 if J2y nnxi^ ~ Vx) > 0, 
fxiv) = otherwise. This has recently received a renewed attention, espe- 
cially since it has been proposed as a model for strong glasses [25] . It is easy 
to check that the presence of a single vacancy in this model allows one to 
empty the whole lattice. 

Among cooperative models we recall FA/ on an hyper-cubic lattice of dimen- 
sion d with 2 < f < d [14]. Here the constraint requires that at least / of the 
surrounding sites are empty in order for the birth/death rate to be non zero. 
As can be directly checked, for all these models it is not possible to devise 
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a finite seed of vacancies wliicli allows emptying the whole lattice. Consider, 
e.g., the case d = 2, f = 2 (with periodic boundary conditions) and focus 
on a configuration which contains two adjacent rows which are completely 
filled. These particles can never be erased, even if the rest of the lattice is 
completely empty: it does not exist a finite defect which can destroy them, 
thus the model is cooperative. The restriction on / comes from the fact that 
all cases f > d are very special because at any value of p there exist finite 
sets of forever blocked particles, namely a fraction of the system is frozen at 
all densities. Although they might be interesting in themselves, they do not 
seem suitable to describe the slow dynamics close to glass-jamming transition 
because they are non-ergodic at any temperature/density. 
A different class of cooperative models is the one we introduced in [15], 
which we call Knight models. On a square lattice kinetic constraints are 
the followings: a move on x can occur only if [the two North-East or the 
two South- West fourth nearest neighbors are empty] and [the two North- 
West or the two South-East forth nearest neighbors are empty]. The defini- 
tion of the North-East, South-East, North- West and South- West neighbors 
and an example of the constraints is given in FigH In formulas fx{v) = 
/,(r/)^^-^^/,(r7)^^-^^ where 

/.(r/)^^-^^ = [(1 - VNElJil - VNE2J + (1 - VSWlJil - VSW2J] 

with NEU = X + 2ei + 62, NE2x = x + d + 2e2, SWl^ = x - 2ei - 62, 
SW2x = x-ei- 2e2, ^^1^ = x + 2ei - 62, SE2x = x + Ci - 2e2, NWU = 
X -2ei + 62 and NW2x = x - ei + 2e2- 

Finally, we recall the one- dimensional East model. In this case the constraint 
requires a vacancy on the right nearest neighbors, that is fx = {1 — Vx+i)- 
Note that on a finite lattice the presence of a single vacancy on the rightmost 
site allows to empty the whole lattice. However the model does not belong 
to the above defined non-cooperative class, since the vacancy should occur 
in a specific position. Indeed, it is usually classified as cooperative: due 
to the direct nature of constraints the relaxation involves the cooperative 
rearrangements of large regions as p ^ 1 [18-20]. 

2.2 Kinetically Constrained Lattice Gases (KCLG) 

KCLG are endowed with a conservative dynamics. A particle in x attempts 
at a fixed rate to jump to a random nearby empty site y and the move 
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a) b) c) 

Figure 1: a) Site x and the four couples of its North-East (NE), South-East 
(SE), North- West (NW), South- West (SW)neighbours. b) Filled (empty) cir- 
cles stand for filled (empty) sites. Here x is occupied and cannot be emptied 
because in the NorthEast-SouthWest direction none of the two couples is 
completely empty, c) Same configuration as in b) except for site x + 2 61+62 
which is empty. Here x can be emptied since both its NW and both its NW 
neighbours are empty, which guarantees 77 e Ax- 

occurs with rate 77^.(1 — riy)fx,y{v)- Here fx,y{v) does not depend on the 
configuration on x and y. Since dynamics preserves the number of particles, 
detailed balance on finite volume is satisfied w.r.t. the measure which is 
uniform on configurations with fixed particle number. On the other hand, on 
infinite volume detailed balance holds w.r.t any Bernoulli product measure 
fip. Again, we can classify KCLG as non-cooperative and cooperative models. 
For the former it is (for the latter it is not) possible to construct a finite group 
of vacancies such that for any configuration it can be moved all over the lattice 
and any jump of a particle to a neighboring empty site can be performed when 
the particle is adjacent to the empty cluster. We will call this mobile cluster 
which facilitates jumps a macrovacancy. As can be immediately checked the 
model in which there are no further kinetic constraints besides hard core, the 
so called normal lattice gas or symmetric simple exclusion process (SSEP), 
is non cooperative and the minimal macrovacancies are simple vacancies. 

Among cooperative models we recall Kob Andersen model (KA) [17] on 
a cubic lattice. Here a particle can jump to a neighboring site only if both 
in the initial and final position at least m of its nearest neighbors are empty, 
where m = 3. Analogously, one can define KA models on hyper-cubic lattices 
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of dimension d and for different values of the parameter m, with 2 < m < d. 
The restrictions on m are due to the fact that m = 1 corresponds to SSEP 
while a model with m > d has a finite fraction of particles that is frozen at 
any density (cf. FA models with / > d + 1). Again, one can directly check 
that all these models are cooperatives according to the definition above. 
Consider, e.g. KA in d = 3 with m = 3 and focus on a configuration 
containing a completely filled slab which spans the lattice and has a two by 
two transverse section. Any finite cluster of vacancies can never overcome this 
filled structure, therefore it is not possible to construct a finite macrovacancy 
that moves everywhere in the lattice. Another possible choice of cooperative 
KCLG are conservative Knight models [15] defined as follows. A move from 
X to y occurs only if the configuration satisfies the requirement needed to 
allow the move in x and in y for the non conservative Knight model defined 
in previous section (plus y should be empty). 

On the other hand, an example of a non cooperative KCLG is provided 
by KA model with m = 2 on a triangular lattice. Indeed, two neighboring 
vacancies form the required mobile macrovacancy as will be further explained 
in Section m 

3 Glass-jamming transition in KCM 

3.1 Ergodicity breaking and glass-jamming transitions 

Let us take the point of view that glass and jamming transitions observed 
in experiments are not just cross-overs but they do correspond to a real 
thermodynamic or dynamic phase transition. For simplicity we will focus on 
lattice models assuming that the presence of the underlying lattice does not 
change qualitatively the physics. 

As explained in the Introduction, from the experiments it is clear that 
ideal glass-jamming transitions (if they exist) have peculiar features com- 
pared to standard first or second order phase transitions. In particular, they 
are characterized by a diverging timescale although no growing static cor- 
relation lengthscale has ever been found, contrary to second order phase 
transition where the relaxation timescale diverges because long-range ther- 
modynamic order sets in. This strongly suggests that ideal glass-jamming 
transitions are purely dynamic phase transitions, i.e. they are characterized 
by an ergodicity breaking without any singularity in the thermodynamics. 
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Is it possible? This question has been addressed in the mathematical physi- 
cal literature (see [42] for a detailed discussion) focusing on interacting hard 
core particle systems on finite dimensional lattices. If particles interact via a 
short-range potential and the rates for exchanging a particle with a nearest 
neighbor vacancy are always strictly positive then it has been proved that in 
two dimensions the ensemble of the Gibbs measures, Q, coincides with the 
ensemble of stationary (under the dynamical evolution) measures. In higher 
dimensions it has been proved that Q coincides with all reversible stationary 
measures, where reversible measure means (roughly speaking) that it verifies 
detailed balance. It is not known if there exist stationary non reversible mea- 
sures. It seems unlikely because they would lead to circulating probability 
currents without any net injection of energy inside the systems, see [42]. In 
any case they are certainly not relevant for the problem of the glass transi- 
tion. Thus, as far as glass-jamming transitions are concerned, one can safely 
assume (as it has been proved in two dimensions) that the ensemble of Gibbs 
measures coincides with the ensemble of stationary measures. Therefore, if a 
dynamical transition takes place, this would lead to more than one stationary 
measure and as a consequence to more than one Gibbs measure: dynamical 
transitions has to be accompanied by thermodynamic transition. 
Would this mean that the ideal glass-jamming transitions must be only of 
thermodynamic origin and no pure dynamical transition can take place? This 
is puzzling, because the absence of any indication of growing static correla- 
tion lengths seems to be at odds with the very existence of a thermodynamic 
transition inducing a diverging relaxation time. If one does not want to 
let down the assumption that experimental glass-jamming transitions cor- 
respond to real phase transitions there remains mainly two possibilities. A 
thermodynamic transition indeed takes place but it is of a completely new 
type and such that is not visible in any given ra— point correlation function. 
A good candidate for that is the Random First Order Theory [7]. Otherwise, 
some hypothesis that lead to the previous conclusion (no purely dynamical 
transition) must be violated. The only one that is reasonable to violate on 
physical grounds is the fact that the rates for exchanging a particle with a 
nearest neighbor vacancy have to be strictly positive. Consider for example 
hard spheres with Brownian or Monte Carlo dynamics (that is a reasonably 
good model for colloidal or even granular systems displaying glass-jamming 
transitions). If the density is high enough some moves can have rate zero. 
This can have a dramatic consequence on the dynamics because now the con- 
figuration space can be broken into subsets that are not connected by any 
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Figure 2: The configuration space of tlie two discs in tlie box is broken 
in disconnected components. It is not possible to reacli dynamically the 
configuration with the disk A in the position of B and viceversa. 

dynamical move when density is large enough, see Fig|21for a trivial example. 
Formally, the Markov chain associated to the dynamical evolution becomes 
reducible contrary to the case of the lattice gas with strictly positive rates. 

As a conclusion the only possibility of having a purely dynamic transition 
seems to boil down to having zero (i.e. degenerate) rates depending on 
the configuration around a site. This encodes the physical effect explained 
previously, see also Fig. 12 that is certainly present for systems of hard 
objects as colloids and granular media. The dynamical transition that may 
be induced by degenerate rates is a reducible-irreducible phase transition. 
At low density (high temperature) almost any equilibrium configuration is 
contained in the same ergodic component. Instead, at high density (low 
temperature), the configuration space sampled with the equilibrium measure 
is fractured in many different components (and no one of them covers almost 
all the space). KCMs are the simplest model that take this phenomenon into 
account. Because of their simplicity they allow to study in great detail the 
issue of glass-jamming transition as purely dynamical transition. 

Indeed for KCMs (see Section 12)) some rates are zero, hence their name 
kinetically constrained, and this implies that the configuration space on fi- 
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nite lattices is generically reducible. In particular, there exist configurations 
which cannot be connected one to the other using allowed moves. For exam- 
ple, for FA model with / = 2 and d = 2 (see Section \2.1\i . a configuration 
which is empty with the exception of two adjacent rows which are completely 
filled can never be connected to any configuration which does not contain this 
slab. This implies that on any finite volume A the process is not ergodic and 
fipS are not the unique invariant measures. The crucial issue is whether er- 
godicity is restored in the thermodynamic limit and whether this depends 
upon the value of p. 

Two final remarks are in order. First, as shown in [36,40], thanks to the 
product form of the equilibrium measure, for KCMs establishing ergodicity at 
density p corresponds to proving that there exists an irreducible set of config- 
urations^ (i.e. configurations connected one to the others by allowed paths) 
which has unit probability w.r.t. pp in the thermodynamic limit. Another 
way to put it is that the only possible way of breaking ergodicity is by the 
reducible-irreducible transition discussed above. The second remark is that 
this type of transition is very different from the critical slowing down related 
to second order phase transitions. In those cases the ergodicity breaking is 
induced by the thermodynamic limit (the number of degrees of freedom goes 
to infinity and so the relaxation time may go to infinity). Instead, for the 
type of dynamical transition advocated above, finite size systems are typi- 
cally non ergodic and the crucial question is whether in the thermodynamic 
limit the ergodicity is restored and up to which density. 

3.2 Dynamical transition in KCMs 

and jamming- boot strap- kcore percolation 

For non-cooperative models it is immediate to conclude that ergodicity holds 
at any density p < 1. Indeed, the probability of finding at least a defect 
or a macrovacancy for an equilibrium configuration goes to one in infinite 
volume. Starting from this macrovacancy, by definition of non-cooperative 
KCSM, one can empty all the lattice. Therefore, two typical equilibrium 
configurations are almost surely connected by a path (with strictly positive 
rates) in the configuration space. The situation is similar for non-cooperative 

^More precisely, we should identify a sequence of sets Tl on 0,1'"^^' such that J-l 
is irreducible w.r.t. the chosen dynamics on with periodic boundary conditions and 
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KCLG. In this case configurations with a macrovacancy can be connected by 
a path which subsequently performs all nearby particle exchanges by previ- 
ously moving the macrovacancy near the interested sites. 
The case of cooperative KCMs is much more involved and indeed it can lead 
to an ergodicity breaking transition. In general, as discussed previously, in 
order to prove ergodicity one should construct an irreducible set and prove 
that it has unit probability w.r.t. /ip in the thermodynamic limit. For KCSMs 
(with periodic boundary conditions) the relevant irreducible set is composed 
by configurations that can be completely emptied by allowed moves, as it can 
be found easily analyzing these systems at low density. This component can 
be identified by the following deterministic procedure: subsequently remove 
all particles for which the constraint is verified until reaching a completely 
empty configuration or one in which there is a backbone of mutually blocked 
particles. Note that this backbone is uniquely determined by the initial con- 
figuration: it does not depend on the chosen order to erase particles. If the 
backbone is empty, we say that the configuration is internally spanned and, 
by definition, it belongs to the relevant irreducible component. Therefore, 
ergodicity breaking takes place if and only if equilibrium configurations con- 
tain an infinite backbone of blocked sites in the thermodynamic limit. Note 
that, as discussed previously, we will consider only choices of the constraints 
such that a blocked cluster has necessarily to be infinite. 

Thus, the problem of the existence of an ergodicity breaking transition 
for cooperative KCSMs can be reformulated as a percolation transition for 
the final configuration of the above cellular automata. For FA models the 
latter coincides with bootstrap percolation procedure [43,44] or k-core perco- 
lation [45] : a particle is removed if it has less than m neighbors. The results 
in [44, 46] establish that bootstrap percolation converges to a completely 
empty lattice for all p < 1, namely no percolation transition occurs. There- 
fore, it is immediate to conclude that FA/ on any hyper-cubic li-dimensional 
lattice does not display an ergodicity breaking transition at any p < 1 for all 
/ and d. On the contrary, for the cellular automata corresponding to Knight 
models, a jamming percolation transition (as we called it) occurs at a critical 
density pc < 1 a.s proved in [15,37]. Therefore, Knight models in infinite 
volume are ergodic for p < pc and non ergodic for p > pc. In Section IHl we 
will explain the mechanism which induces the percolation transition for the 
blocked structures of Knight models and discuss its character. 
For cooperative KCLG, the prove of ergodicity is more involved. For exam- 
ple, for KA model on square lattices with m = 2, we have shown that the 
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irreducible component which has unit probabihty in the thermodynamic hmit 
is the one composed by configurations which can be connected by an allowed 
path to a configuration which have a frame of vacancies on the boundary [36]. 
Establishing that this set has unit probability in the thermodynamic limit 
is more involved than in the corresponding KCSM (i.e. FA model). This is 
due to the fact that there do not exist a deterministic bootstrap-like proce- 
dure which allows to establish whether a configuration belongs or not to the 
irreducible set. 

Finally note that, apart from ergodicity breaking, other types of transi- 
tion might take place in lattice models for glass-jamming transitions: (1) a 
change from an exponential to a stretched exponential relaxation for density- 
density correlation in KCSM or its Fourier transform for KCLG, (2) a dif- 
fusive/subdiffusive transition in the behavior of a probe particle and (3) a 
break down of the conventional hydrodynamic limit. 

In the case of cooperative KCSM numerical simulations [17,29-31] seems 
to indicate that connected correlation functions, e.g. < 'r]x{t)'r]x{^) >p — < 
Vxit) >p< Vxi^) >p, decay in time as stretched exponential, exp[— (t/r)^], 
with a stretched exponent f5 that is less than one and density dependent 
{P decreases as p is increased). Analytical works based on mode coupling 
approximations [28] support these results. The standard explanation of this 
behavior is dynamic heterogeneity, which has indeed been experimentally de- 
tected [4] and consists in having different relaxation times in different regions 
of the system. The superposition of these exponential relaxations would then 
lead to an effective stretched exponential behavior. 

Although this is likely what happens for KCSMs and real liquids, it has 
been recently proved in [23] that relaxation is exponential at very large times 
for all cooperative and non-cooperative models (in the ergodic regime) con- 
sidered so far. More precisely, relaxation is exponential on times of the order 
of the inverse of the spectral gap of the Liouvillian operator generating the 
dynamics (that is the slowest relaxation time over all one time quantities) 
and this is proved to be finite at any p in the ergodic regime. The dis- 
crepancy between rigorous and numerical results is likely due to the different 
time regimes that are investigated: a correlation function can be well approx- 
imated as a stretched exponential for times t of the order of its relaxation 
time ^ r and decay as an exponential when t » t. An interesting question 

■^Here we mean the time over which the specific normalized connected correlation func- 
tion has decayed to e^^, which may of course be shorter than the slowest relaxation time 
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is whether the effective numerical behavior coming from fit on times of the 
order of r can be made precise: in the scahng hmit t oo, with t/r of the 
order of one, do correlation functions converge to stretched exponentials?. 

Another transition different from the ergodicity breaking is a diffusive/ 
subdiffusive transition for the motion of tagged particles. The result that 
the self-diffusion coefficient, -Ds(p), is positive at any density less than one, 
which is known to hold for SSEP cannot be immediately extended to KCLG 
due to the degeneracy of jump rates. In particular for KA models, from 
numerical simulations [17] it had been conjectured that Ds > only below 
a finite critical density. However, in [40] it has been proved Dg > a.t any 
density p < 1 for non-cooperative models and the result have been extended 
to cooperative models (in the ergodic regime) in [41]. 

Finally, a third possibility is a breakdown of the conventional hydrody- 
namic limit that usually holds on long length and time scales [42]. For KCLG, 
if the constraint is released, the macroscopic density profile evolves via the 
diffusion equation dtp = V(-D(p)Vp) where D{p) = 1 for the normal lattice 
gas. Again, due to the presence of constraints, one cannot apply the tech- 
niques which have been developed to establish this hydrodynamic limit for 
stochastic lattice gases [47]. Furthermore, for some cooperative KCLG, it has 
been conjectured [50] that the macroscopic diffusion coefficient D{p) would 
vanish at high density, leading to a sub-diffusive evolution of density profiles. 
In [48] it is proved that for a class of non cooperative models hydrodynamic 
limit holds with D{p) vanishing (as power law of 1 — p) only for p — > 1. 
Extending these results to cooperative models would require establishing the 
scaling of relaxation times on finite lattices with the size of the lattice, which 
has not yet been settled for cooperative KCLG. 

4 Non-cooperative models as renormalized lat- 
tice gases 

In Section 121 we explained that, thanks to the presence of the macrovacan- 
cies (for KCLG) and defects (for KCSM), in the thermodynamic limit non 
cooperative models are ergodic at any p < 1 as it occurs for the correspond- 
ing models without kinetic constraints. In this Section we shall explain that 
the slow dynamics for non cooperative models can be understood in terms 

given by the inverse of the spectral gap. 



14 



Figure 3: The triangular lattice. For KA with m = 2 a particle in x can be 
moved to the neighbouring empty site x + (ei + e2)/2 thanks to the fact that 
the other empty site {x + (— ei + e2)/2) is nearest neighbors of the particle 
both in its initial and final position. 

of motion of macrovacancies or defects, i.e. as renormalized lattice gases. 
Indeed these mobile regions, when p — > 1, essentially perform independent 
random walks and substantial relaxation takes place only when they pass by. 
So the high density dynamical behavior is encoded in their properties (size, 
density, timescale for motion, . . . ). The analysis and the physical scenario 
is very similar for non-cooperative KCSM and KCLG. In the following we 
shall focus on KCLG which are in a sense richer because, apart from the re- 
lation timescale, one can study directly the self-diffusion coefficient of tagged 
particles, i3's,bulk diffusion coefficients, -D(p), etc. 

Let us first focus on the scaling with density of Dg. Consider for example 
KA model with m = 2 on a triangular lattice (see Section 12)). The triangular 
lattice A, represented in Fig. is the union of sites in a square lattice 
Ai and in its dual A2 which is obtained by displacing Ai of (ei -|- e2)/2 
with ei = (1,0), 62 = (0,1). Two sites {x,y} G A are nearest neighbors if 
X — y = ±ei or X — y = ±(ei + e2)/2 or x — y = ±(— ei + 62) /2. Dynamics 
can be reformulated in terms of vacancy motion: a vacancy can move to a 
neighboring occupied site only if it has at least one empty neighbor both in 
its initial and final position. Since two neighboring sites of the triangular 
lattice share a common third neighbor, it is immediate to see that any of 
the two vacancies of a neighboring couple can move to the common third 
neighbor. Therefore a couple of neighboring vacancies is a finite size freely 
mobile cluster. Furthermore it is possible to move any given particle into 
an empty nearest neighbor provided the couple of vacancies is nearby. For 
example, if we want to move a particle from x to x + (ei +e2)/2 it is sufficient 
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to put the macrovacancy in x+(ei+e2)/2, ei+e2)/2, as shown in FigOl 

Thus, a couple of neighboring vacancies is a macrovacancy according to the 
definition in Section|21 A simple heuristic argument based on the independent 
motion of these macrovacancies leads to the correct high density dependence 
of Dg. Let us focus on dimensions larger than two (for lower dimension the 
reasoning changes due to recurrent properties of random walks) . Call pd the 
density of macrovacancies and the timescale on which they move. The 
self-diffusion coefficient is expected to be proportional to the inverse of the 
time Tp on which each particle moves of one step. On the timescale the 
number of particles that have jumped is of the order V pd where V is the total 
number of sites. Thus we find 

^Fprf oc Vp. (1) 

Td 

As a consequence, at density close to one, we get Ds oc Pd/Td- Since for 
the m = 2 KA on a triangular lattice macrovacancies are formed by two 
neighboring vacancies, in the limit p ^ 1 we get pd oc (1 — p)^ and Tfi oc 0(1). 
Hence, Ds oc (1 — p)^. The same arguments leads to Ds oc (1 — p) for 
SSEP, where macrovacancies are single vacancies. The result for SSEP has 
been rigorously proved in [49] by establishing upper and lower bounds. For 
KA model on a triangular lattice, as we have shown in [41] (see also [40] 
for a different choice of non cooperative constraints), it is also possible to 
turn the heuristic argument into a proof deriving upper and lower bounds 
q(1 — p)^ < -Ds(p) < c„(l — p)^ with Cu and q independent from p. 

The representation of non-cooperative KCM as renormalized lattice gases 
turns out to be useful also to obtain these rigorous results. Let us recall the 
basic strategy since it allows to further understand the role of macrovacancies. 
The idea is to start from the proof for SSEP in [49] and to modify it by letting 
macrovacancies play the role of vacancies. The procedure is the following: 
first construct a suitable auxiliary model that is easy to analyze and then, 
using a variational formula for Z)s, show that Ds > h{p)D^^^ where h{p) is 
a function determined explicitly (that is strictly positive for p < 1). The 
technique for the upper bound is more classical [49] and the lower bound 
is more important because it shows that, despite the dynamics slows down, 
particles still diffuse on large distances and timescales. Therefore, we will 
just focus on the lower bound in the following. 

In [49] the auxiliary model is constructed in this way: the tracer has a 
neighboring vacancy at time zero and the only possible moves are jumps of 
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the tracer and exchanges of the occupation variables in y and w, where both 
y and w are nearest neighbors of the tracer. Note that the latter moves are 
not allowed for SSEP but one can reconstruct them trough a path of allowed 
moves thanks to the presence of the vacancy. For the auxiliary model it is 
immediate to prove that D'^'^^ > and is bounded from above and below 
by density independent constants. Indeed, a move of the tracer from x to 
X + Ci can always occur via two steps: since the tracer has always at least 
one neighboring vacancy, this can be brought (in one move) in x + ei and 
then the jump of the tracer from x to x + ei can occur. Using a variational 
formula for Ds [49] and comparing the dynamics of the normal lattice gas 
and the auxiliary model leads to D^^^^ > c (1 — p) Df^^. The term (1 — p) 
comes from the requirement of having at least a vacancy at time zero near the 
tracer; c is a density independent constant which accounts for the maximal 
length of the path needed to exchange the occupation variable in x and w 
trough neighboring jumps (the only moves allowed for SSEP). For KA on the 
triangular lattice one defines the following auxiliary process. We require a 
macrovancancy near the tracer in the initial configuration and we allow only 
moves which either displace the tracer to a nearest neighbor or move the 
macrovacancy into another couple of sites which are also near the tracer (for 
a precise definition of the auxiliary model see [41]). Again, the latter move 
is not allowed in KA but can be reconstructed by a finite path of allowed 
moves. As before, it is easy to show that _D^"^ > and is independent from 
p. Comparing KA and auxiliary dynamics we get Dg > c (1 — p)^, where 
(1— p)^ comes from the requirement of having a macrovacancy. Note that this 
procedure is generalizable to all non cooperative models: if the macrovacancy 
is formed by q empty sites it leads to Ds{p) > c(l — pY- Therefore, at least 
for the motion of the tracer, non cooperative KCLG are simply renormalized 
SSEP: the typical diffusion time of the tagged particle (i.e. the inverse of 
Ds) goes here as the inverse of the density of macrovacancies, instead of 
vacancies. 

The strategy of comparing dynamics with the one of a faster uncon- 
strained model and reconstructing the moves of the latter by a proper path 
of allowed moves can be also used to derive upper bounds for the density and 
size dependence of the relaxation time, r, for non-cooperative KCLG and 
KCSM (see [40] and [23], respectively). The paths are always constructed 
using the fact that, if there is a macrovancancy (defect), we can move it ev- 
erywhere and facilitate any nearest neighbor jump (birth/death) of particles. 
Therefore comparison with the unconstrained model in general gives an extra 
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factor (1 — pY (the cost of creating the mobile region) times a term related 
to the length of the path (which can bring a dependence both on the lattice 
size, L, and on p). Note that in principle also an entropy term should be 
accounted (for the bounds on r and Dg), counting all the possible configu- 
rations which have to pass trough the same bottleneck [23,36,40]. However, 
this is just a constant factor for non cooperative models (more precisely it 
is the entropy on a number of sites proportional to the size of the macrova- 
cancy) and cannot change the scaling in 1 — p. 

From results in [23,40], one has that for non-cooperative models the depen- 
dence of r on the lattice size L is the same as for models without constraints: 
T oc L'^ for KCLG, r bounded by a finite constant uniformly in L for KCSM 
(where r is the inverse of the spectral gap of the LioviuUian operator, see [42]). 
On the other hand, since relaxation occurs trough the macrovacancies which 
become rarer at higher density, the density dependence is different. For exam- 
ple, for KA on a triangular lattice with particle sources at the boundary, we 
find L^/(l — p)q < t < L^/(l— p)^Cu, to be compared with ci L'^ < t < Cu L'^ 
for SSEP. 

5 Jamming percolation on Bethe lattices 

Bethe lattices (or in the mathematical literature random c-regular graphs) 
are often used in the physics literature as an approximation of finite dimen- 
sional, e.g. hyper-cubic, lattices. Because locally they have a Cayley tree 
structure with connectivity k + 1, their analysis is greatly simplified and one 
can often obtain exact results on the thermodynamics and the dynamics of 
the model embedded on them. For proper choices of the constraints, KCMs 
display a jamming transition on Bethe lattice [35,51] and, thus, provide an 
almost solvable example of jamming transition. In the following we present 
a summary of the results that have been obtained focusing on a simple case: 
the f-facilitated Fredrickson- Andersen model on a Bethe lattice with connec- 
tivity k (we consider A; > / > 1 in order to avoid any finite blocked cluster). 
The results remain qualitatively the same for particle models but they are 
more tedious to derive [36]. 

As discussed previously, a jamming (and ergodic) transition for FA models 
will take place if and only if at a certain density an infinite cluster of blocked 
(jammed) particles appears. All the particles inside this cluster must have 
at least / blocked neighbors in order to be blocked themselves. Thus, this 
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Figure 4: Solution of the equation Q as a function of p for /c = 3 and / = 2. 

infinite cluster is the spanning cluster of the so called /—core percolation [45], 
(or bootstrap percolation [52]) a problem that has received a lot of attention 
recently [53-55]. 

For a given site the Bethe lattice local structure is Cayley-tree like with 
k branches going up from each node and one going down. Using this crucial 
feature it is easy to write a self-consistent equation on the probability P that 
a site is occupied by a blocked particle (belongs to the /—core) because it 
has more than k — f + 1 neighboring particles above it which are blocked 
(without taking advantage of vacancies below) [52]: 



Except ioT k = f for which the transition is continuous, the solution of this 
equation leads to P = for p < pc and P = Pc + 0{^ p — pc) with Pc > 0, see 
the example for k = 2> and / = 2 in Fig. |3] The number of blocked particles 
or the number of sites belonging to the /—core are polynomial functions of 
P that can be derived easily and that have the same behavior of P. They 
are at the same time discontinuous as in first order phase transitions and 
singular because of the square root, as in second order phase transitions. 
The mechanism behind this behavior has been understood in detail [53, 
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54]. In particular the singular square root behavior is due to the extreme 
fragility of the infinite spanning jammed cluster at the transition and to the 
existence of a related diverging lengthscale. In fact, close to the transition, a 
given jammed particle is connected to clusters of jammed particles that are 
only marginally blocked, i.e. such that removing only one of their blocked 
neighbors is enough to unblock them. These clusters form the so called 
corona of the infinite jammed cluster (or /—core spanning cluster). The size 
distribution of corona clusters can be computed analytically [53] and it has 
been found that their average size diverges at the jamming (/—core percola- 
tion) transition as l/-\/p — pc coming from the jammed phase. This explains 
the singular square root behavior found previously [53]: roughly speaking, 
decreasing the temperature from Pc + e to, say, pc + e/2 (where < e << 1) 
one unblocks first a number of particles proportional to Ne. However, un- 
blocking a particle unblocks also all the marginally blocked (corona) particles 
attached to it. This leads to a "domino effect" such that the effective number 
of unblocked particles is A^e x l/-\/e, i.e. the net change in the fraction of 
blocked site is y/ p — pc- Since corona clusters, i.e. marginally stable clusters, 
percolate at the transition one expects an associated diverging lengthscale. 
This is indeed the case, see [53,54]: coming from the jammed phase there is 
lengthscale that diverges as (p — pc)~^^^. 

The analysis of the jamming transition coming from the unjammed phase 
is much more difficult and one has to resort to numerical simulations. In 
[34,51] the persistence, 

i 

where Vti{t) equals one if the site i has remained in the same state (empty 
or occupied) from time to time t and zero otherwise, and the occupation 
variable correlation function, 

i 

have been measured starting from equilibrium initial conditions. We plot in 
Fig. Elthe persistence as a function of time for the f = 2, k = 3 case. Different 
curves corresponds to different densities approaching the critical point ^. The 

■^Note that in the work [51] we used a different terminology: spin instead of occupa- 
tion variables and temperature instead of density. The relationship between density and 
temperature T is p = 1/(1 + e^^/^). 
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Figure 5: Persistence as a function of time for the f = 2, k = 3 case for 
different temperatures approaching the critical temperature ~ 0.48. 



correlation has a similar behavior [34] . The existence of the plateau is a direct 
consequence of the discontinuous character of the jamming transition: the 
fraction of blocked particles is strictly positive at the transition coming from 
the jammed phase instead it is zero coming from the unjammed one. This 
translates into a dynamical behavior such that Ymvpyp^ limj_^oo P{t) = and 
limt^oo ^^"^pypc P{^) > and equal to the fraction of blocked sites. Associated 
to the non-commutation of these two limits there is a diverging time scale, 
the time on which P{t),C{t) equals, say, half of their plateau value. The 
numerics indicates that it follows a power law divergence r oc [pc — p)~^ with 
7 ^ 2.9. 

Finally, it has been found a diverging dynamical lengthscale associated 
with the divergence of the timescale [51] (let us recall that no static correla- 
tion can be found because the equilibrium measure is uncorrelated from site 
to site). In order to unveil this length one has to focus on the fluctuations of 
the persistence or the correlation [57]. These are encoded in the dynamical 
susceptibility 



X{t) = N{ 



)/T 



x(t) develops a peak (see FiglHl ) for t oc r that diverges for T Tc. This 
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Figure 6: Equilibrium dynamical susceptibility x(t) vs time t at tempera- 
ture T. System size = 2^'^. 

divergence is related to the fact that close to Tc more and more particles have 
to evolve in a correlated way in order to make the system relax. Note that, 
exactly as for ferromagnets the divergence of the fluctuations of the magne- 
tization is a signature of a diverging lengthscale, so it is the divergence of the 
peak of x{t) at the jamming transition. However, the precise characteriza- 
tion of this dynamic lengthscale and its relationship with the results on the 
jammed phase have still to be worked out. 

Let us summarize the features of the jamming transition found on the 
Bethe lattice. It has a first order character because the fraction of jammed 
particles and the long-time limit of correlation and persistence functions 
(called in spin-glasses the Edwards-Anderson parameter) are discontinuous 
at the transition. At the same time it is characterized, as second order phase 
transitions, by a diverging timescale and diverging lengthscales. Although 
the detailed connection between time and lengthscales (from the jammed and 
unjammed phase) has still to be worked out from numerics it seems clear that 
all these quantities diverge as power law in \p — pd- It is interesting to note 
that qualitatively these features are also shared by mean-field disordered sys- 
tems, as p-spins models, close to their dynamical transition and the Mode 
Coupling Theory of the glass transition [6]. The following section is devoted 
to the analysis of jamming transitions in finite dimensional models. As we 
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will see some of the features found on the Bethe lattice persist, for example 
the mixed (first and second order) character of the transition. Instead oth- 
ers, as the divergence of the relaxation times, change and actually become 
closer to what expected for glass-jamming transitions: now the logarithm of 
r diverges as a power law in |p — Pc|, i-e. in a Vogel-Fulcher like form. 

6 Jamming percolation on finite dimensional 
lattices 

In this Section we give an extended explanation of the results obtained in 
collaboration with D.S.Fisher [15,37] for the Knight models defined in Section 
121 As already mentioned an ergodicity breaking transition occurs at pc < 1: 
above pc a finite fraction of the system is frozen. This is due to an underlying 
percolation transition for the clusters of mutually blocked particles, which we 
called jamming percolation. 

Let us give a brief summary of the main results letting for the next sections 
their derivation. 

As discussed in Section 3 the dynamical transition of the Knight mod- 
els can be studied using the following cellular automata [37]. Start from an 
initial configuration sampled with equilibrium measure, /ip, then empty an 
occupied site if the constraint is satisfied. Then continue the procedure until 
no more particles can be removed. This final configuration is either com- 
pletely empty or one that contains a percolating cluster of particles which do 
not satisfy the constraint. Let us call poo the density of this final configura- 
tion. We will show that at the critical density of site directed percolation on 
a square lattice, pdp — 0.705, a discontinuous percolation transition occurs: 
Poo = for p < and poo > for p > p^p. Since the final backbone for 
this cellular automata contains all the particles that are frozen under the 
stochastic evolution of Knights, we conclude that below p^p Knight model 
is ergodic and ergodicity is broken above. Furthermore, the fraction of the 
system which is frozen coincides with p^ and has a finite jump at the tran- 
sition. 

As for the FA on a Bethe lattice, we find that the dynamical transition tran- 
sition has a first order character with a discontinuous Edwards-Anderson 
parameter and, at the same time, is characterized by diverging time and 
lengthscales, as second order like transitions. 
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The diverging lengthscale can be determined studying finite size effects. Con- 
sider the model on a finite lattice with periodic boundary conditions and 
evaluate the probability R{L, p) that it can be completely emptied by allowed 
moves, i.e. the probability that is internally spanned. Since at any fixed 
density p < pdp the final backbone of the cellular automata is empty, R{L, p) 
goes to one in the thermodynamic limit. On the other hand, if the limit 
P y Pdp is taken first, R{L,p) goes to zero. The relevant length H(p) is the 
crossover length which separates the two regimes: \imL^oo,p/'pd_p R{L, p) = 
for L/E{p) — > and \imL^oo,p/'pdp p) = 1 for L/E{p) oo. In other 
words, is the dynamical length which corresponds to the typical size 
of blocked clusters (i.e. the incipient percolating clusters for jamming per- 
colation). As we will explain, by analytical arguments we find [15] that 
log S(p) ^ (p — pdp)~^- Here p = z/||(l — z) where z/y ~ 1.73 is the crit- 
ical exponent of the parallel directed percolation (DP) correlation length, 
,^11 ^ (p — pdp)^'^", and z ~ 0.63 is the exponent relating parallel and trans- 
verse DP correlation lengths, ~ We also rigorously proved ^ the upper 
and lower bounds Ci exp(c2^||~^) < S(p) < C3exp(c4^j|) [37]. These estabhsh 
that S diverges faster than power law and the lower bound coincides with 
the value we expect for H. 

Due to this diverging length, the relaxation time r should diverge as with 
z >2. Indeed, relaxation should require the diffusion of regions of density S 
(see [36] for similar results for FA model). In [23] the rigorous bound r > H is 
proved, where r is the inverse of the spectral gap of the LioviuUian operator C 
(i.e. the worst relaxation time on all one time quantities). Therefore Knight 
models exhibit a Vogel-Fulcher like relaxation: the log of times diverge as 
power law when the critical density is approached. This is different from 
the findings for FA model on Bethe lattices, were r diverges as power law 
(see Section EJ. Since in the Knights model there are no static interactions 
between the particles it is clear that their glass-jamming transition is purely 
dynamical and not related to any thermodynamic transition. 

In the following we shall explain the arguments leading to the proof of an 
ergodicity breaking transition at pdp and of its mixed first / second order char- 

'*The existence of the two different correlation lengths, which are due to the asymmetry 
of DP and the power law divergence of , are given for granted in physical literature. In- 
deed, they have been verified both by numerical simulations and analytical works trough 
renormalization technique (see [56] for a review). However, a rigorous mathematical proof 
of these results for DP is still lacking. Our bounds for S are rigorous modulo the assump- 
tion ix_ = Cjr, z < 1. 
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Figure 7: a) Sites inside circles are those obtained starting form the origin 
and choosing its NE and SW neighbours and the NE and SW neighbours of 
the latter and so on. b) Sublattice obtained erasing all sites from the square 
lattice except those inside circles c) Oriented sublattice: arrows go from each 
site to its two NE neighbours (w.r.t. the original lattice) 

acter (we refer to [37] for rigorous proofs). These will unveil the mechanism 
which induces these remarkable properties which are typical of an ideal glass 
transition. Also, it will be clear why we choose the specific form of Knight 
constraints and how one can modify or generalize them to higher dimension 
keeping a transition with the same character. 

6.1 Ergodicity breaking for p > p^p 

Let us prove that ergodicity is broken above pdp, namely a finite fraction of 
sites is frozen. This is equivalent to showing that the origin belongs with 
finite probability to a percolating cluster of frozen sites. Here, as in [15], 
we call a site frozen if it cannot be unblocked even by first emptying with 
allowed moves an arbitrarily large number of sites. Consider the sublattice 
which is obtained from the square lattice by erasing all sites except the origin, 
its North-East and South- West neighbors, the North-East and South- West 
neighbors of the latter and so on as in Fig. [7|a), b). Then construct a directed 
graph on this sublattice by drawing arrows connecting each site to its North- 
East neighbors as in Fig. He). Notice that this corresponds simply to a two 
dimensional square lattice (only rotated and squeezed). Also, since there are 
no static correlations in the equilibrium measure, if the configuration on the 
original lattice is chosen with fip, the same holds for the configuration on the 
sublattice. Therefore the results for site Directed Percolation [56] imply that. 
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if p > pdp, the origin belongs with finite probabihty to an infinite directed 
(i.e. following the direction of the arrows) percolating cluster of occupied 
sites. It is now easy to check that if we restore the whole lattice and consider 
Knight dynamics, all sites belonging to the cluster which percolates on the 
sublattice are frozen: a finite fraction of the system is frozen above pdp- 
Indeed, each of these sites has at least one occupied NE neighbor and at least 
one occupied SW neighbour, therefore it is blocked along NE-SW diagonal. 
In the following, we call NE-SW cluster a directed occupied path on the 
sublattice obtained with the procedure above. The definition is analogous 
for the NW-SE cluster, where the sublattice is the one constructed erasing 
all sites but the origin, its NW and SE neighbours and so on. Note that 
in order to establish the existence of frozen clusters we have not used the 
possibility of blocking along the NW-SE diagonal. On the other hand, this 
should be taken into account if we want to show that a given configuration 
does not contain frozen sites as we will do in the following section. 

6.2 Ergodicity for p < p^p 

Let us prove that the system is ergodic in the thermodynamic limit for p < 
Pdp, which corresponds to showing that the fraction of frozen sites is zero. 
Note that if the rule contained blocking only along one of the two diagonals, 
ergodicity would follow immediately from the fact that occupied directed 
paths do no percolate below pdp- However, for the constraints we have chosen, 
a percolating directed path implies a frozen cluster but the converse is not 
true. Indeed, due to the fact that a site can be blocked along either the 
NE-SW or the NW-SE diagonal or both, a NE-SW cluster can be blocked 
either if it spans the lattice or if it is finite but both its ends are blocked by a 
T-junction with NW-SE percolating paths (see Fig. |Hla)). By using such T- 
junctions it is also possible to construct frozen clusters which do not contain 
neither a spanning NE-SW nor a spanning NW-SE cluster: all N15-SW (and 
NW-SE) clusters are finite and are blocked at both ends by T junctions by 
finite NW-SE (NE-SW) ones (see figure |Hlb)). As we will show in detail, 
these T-junctions are crucial to make the behavior of jamming percolation 
very different from site directed percolation although the two transitions 
share the same critical density. 

Let us now show how the proof or ergodicity below pdp works. The strat- 
egy consists in constructing a set of configurations, J-'l, on and show that 
J^L are internally spanned and lim^^^oo /^p(-^l) = 1- This, as explained in 
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Figure 8: a) A NE-SW non spanning cluster blocked by T-junctions with 
two NW-SE spanning clusters; b) Segments parallel to the NE-SW (NW- 
SE) diagonal stand for NE-SW (NW-SE) clusters. The depicted structure is 
blocked, even if it does not contain neither a NE-SW nor a NW-SE spanning 
cluster 

Section |S1 concludes the proof of ergodicity for p < pdp and (together with 
results in previous section) establish that ergodicity breaking occurs at pdp. 

Consider a configuration within which there is an empty square of size 
n X n and focus on the sufficient conditions to empty the next shell, i.e. to 
construct an allowed path which empties the n + 2 x n + 2 square. The initial 
vacancies guarantee that we can empty a centered segment of length n — 4 
external to each side. Consider for example site x in fig. IHla), both its SW 
and SE couples are inside the empty square: the constraint is satisfied on 
both diagonals and x can be emptied. It is immediate to check that, even 
if the rest of the lattice is occupied, this procedure can be continued until 
emptying the triangular structure above each side depicted in Fig. Elb). 
This is formed by centered lines whose length decreases of four sites at each 
step, therefore its overall height is I{n/4) where I{x) stands for the integer 
part of X. In order to empty the n + 2 x n + 2 square, five sites remain to 
be emptied on each corner: those indicated by question marks in FigElb). 
Consider for example one of these site on the top left corner, indicated as y 
in FigJini The constraint along the NW-SE diagonal is verified, since both 
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Figure 9: a) The empty nxn square. Site x can be emptied since its SE and 
SW neighbours are all inside the empty square, b) The triangular structure 
of height /(n/4) which can be emptied above each side. ? denotes the 5 
sites on each corner which are not guaranteed to be emptied (unless a proper 
condition is fulfilled by the configuration outside the empty nxn nucleus). 



its SE neighbors are inside the empty square. Instead, since neither the NE 
nor the SW couples are contained in the empty region, y can be blocked 
along NE-SW diagonal and in this case we cannot empty it. However, it is 
possible to derive a necessary condition for y to be frozen, which gives in turn 
a sufficient condition to construct an allowed (i.e. verifying the constraint at 
each elementary move) path in the configuration space which allows to empty 
y. Focus on the square of size n/2 + I{n/4) centered on the top left corner 
and containing half of the triangular structure associated to the up and left 
sides, i.e. the square inside the dashed region in Fig. Mp)- Consider the 
sublattice constructed from y inside this square by erasing all sites except its 
North-East and South- West neighbors and so on (as we did for the origin in 
figure [7j). It can be directly checked that y can be frozen only if it belongs to 
a NE-SW path, i.e. a directed path on this sublattice, which spans the square 
(see Fig. fTTH) . This is due to the fact that any occupied cluster along the NW- 
SE diagonal can be unfrozen if it terminates inside the empty region [nxn 
square plus triangular region). Therefore, a sufficient condition to construct 
an allowed path to empty y (and the other remaining four sites at the top left 
corner) is that in the dashed square region there is not a percolating NE-SW 
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Figure 10: Here we depict the top left square in Fig. El b) and show the 
necessary condition for y to be frozen: it should belong to a NE-SW cluster 
spanning the dashed square. Otherwise: if one of its end is free it can be 
unblocked from there; if they are both blocked by T-junctions with NW-SE 
clusters these terminate inside the empty region and can be unblocked from 
there. 
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cluster. Since the size of this region is proportional to n and p < p^p, such a 
cluster occurs with a very small probability, exp(— n/^y) for large n with 
the parallel correlation length for DP [56]. As a consequence the probability 
of emptying y increases at least exponentially fast to one when the size of 
the empty nucleus, n, is increased. The same argument can be applied to the 
other four sites indicated by question marks and to those on the other three 
corners of the lattice, concluding that the probability P{n —>■ n+2) to go from 
the empty n x n to the empty n + 2 x n + 2 nucleus is bounded from below 
by (1 — 20 exp(— n/^||)). Thanks to the exponential increase towards one, 
the probability that this procedure can be continued up to infinite size stays 
strictly positive because it is bounded from below by Yliln -^(^ — ^ i + 2) > 0. 
Finally, using the 0{L'^) different positions for the initial empty nucleus, we 
conclude that the probability for Al to be internally spanned converges to 
one in the thermodynamic limit at any p < pdp. This is due to the fact that, 
in order to prevent the expansion of a large empty nucleus, we should require 
long DP clusters, which are highly improbable below p^p. 

6.3 Discontinuity of the transition 

Results in Sections 16. II and 16.21 prove that an ergodicity breaking transition 
occurs for Knights at pdp-, due to the percolation transition of blocked struc- 
tures. In the present and the following section we show that this transition, 
which we call jamming percolation, have features which are qualitatively dif- 
ferent from those of DP and any conventional percolation transition. Indeed 
it is discontinuous, i.e. the density of the frozen clusters have a finite jump 
at the transition (the critical clusters are compact rather than fractal) and 
their typical size increases faster than any power law when p p^p- 

In order to prove discontinuity we construct a set of configurations which 
have finite probability on the infinite lattice at pdp and show that the ori- 
gin is frozen for all these configurations. We will make use of the blocked 
structures containing T-junctions among NE-SW and NW-SE clusters which 
have been introduced in Section Consider a configuration in which the 
origin belongs to a NE-SW path of length £o/2: this occurs with probability 
Po- Now focus on the infinite sequence of pairs of rectangles of increasing 
size li X lijVl with i\ = ^o? = 2£j_2 and intersecting as in FiglTTl If 
each of these rectangles with long side along the NE-SW (NW-SE) diagonal 
contains a NE-SW (NW-SE) percolating path (dotted hues in Fig. ITD b)). 
the origin is frozen. This can be directly checked: an infinite backbone of 
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Figure 11: a) The sequence of intersecting rectangles of increasing size ii x 
l/12£j describer in the text, b) Dotted non straight hne stand for NE-SW 
(NW-SE) clusters spanning the rectangles with long side in te NE-SW (NW- 
SE) direction, c) Frozen structure containing the origin 

mutually blocked particles constituted by pieces of these paths connected by 
T-junctions (continuous line in Fig. ^Jc)) occurs. Therefore the probability 
that the origin is frozen, g(p), is bounded from below by 

q{p) > Po n ^(^^)' 

i=l,oo 

where P{i-i) is the probability that a rectangle of size li x l/12£j with short 
side in the transverse direction is spanned by a DP cluster. This probability 
converges to one exponentially fast in £j due to the anisotropy of critical clus- 
ters for directed percolation. Recall that there are a parallel and a transverse 
length for DP with different exponents, i.e. a cluster of parallel length I has 
typically transverse length Let us divide the ii x l/12£j rectangle into ^\~^ 
slices of size x 1/12£|. For each slice the probability of having a DP cluster 
along the parallel direction at pdp is order unity. Thus, the probability of not 
having a DP cluster in each of the slice is 1 — -P(^j) = 0[exp(— c£^~^)]. From 
this result and above inequality for g, it follows immediately q{pdp) > 0. 

Therefore the infinite cluster of jamming percolation is "compact", i.e. of 
dimension d at the transition. Note that to obtain discontinuity two ingredi- 
ents of the constraints were crucial: the existence of two transverse blocking 
directions each with an underlying percolation transition; the anisotropy of 
these transitions. Indeed, anisotropy is necessary to have the increase to- 
wards one of the probability that above rectangles are spanned when their 
size is increased. In turn, this is necessary to get a finite probability for the 
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construction which freezes the origin. 

6.4 Dynamical correlation length 

In this Section we explain the arguments in [15] leading to the result log S(p) ^ 
k{p — Pdp)"'^ for p < pdpi where p = z/||(l — 2) ~ 0.64. As already explained, 
S(p) is the crossover length dividing the regime in which the probability for 
a finite region to be internally spanned goes to zero or to one when p pdp, 
i.e. the length below which finite size effects are important. For jamming 
percolation it corresponds to the diverging size of the incipient spanning clus- 
ters and for Knight models to the typical size of the region that has to be 
rearranged to unblock a given site. Let us sketch separately the arguments 
leading to logS(p) > ki{p-pdp)~^ and to logH(p) < ku{p-pdp)~^ with ki, ku 
two positive constants. 

To establish the lower bound we construct a set of configurations contain- 
ing a frozen backbone and we show that the probability of this set goes to 
one when p y pdp and L — >• cx) with logL < ki{p — Pdp)~^- Again, we make 
use of the T-junctions described in Section 16.21 Consider the set of NE-SW 
and NW-SE paths of length s intersecting as in Fig. ^| 

As can be directly checked, this structure can be emptied only starting 
from its border since each finite directed path terminates on T-junctions with 
paths in the opposite direction. Therefore, if the construction is continued 
up to the border and we consider periodic boundary conditions, all sites 
in the structure are frozen. Furthermore, a similar frozen backbone exists 
also if one or more of these paths is displaced inside an adjacent rectangular 
region of size s x s/6, as shown in Fig. ^| Therefore the probability that 
the regions is not internally spanned, 1 — R{L,p), is bounded from below 
by the probability that each of these rectangles contains at least one path 
connecting its short sides. This leads to 

R{L, p) < n{L, s) exp(-cs^-^) (3) 

where n{L,s) = 0{L/sY is the number of the rectangles contained in the 
structure and exp(— cs^~^) is the probability of not having a DP path in 
a region s x s/6 as long as s < (see previous section). The inequality 
(0) leads to Yim.L-^oo,p/pcipR{L^ p) = for exp(c.^y~^)/L 00, therefore 
logS >ki{p- pdp)~^- 

On the other hand, in order to establish the upper bound, we show that 
for logL > ku{p — Pdp)~^ there is typically an nucleus of vacancies which 
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Figure 12: The frozen structure described in the text: continuous hues stand 
for occupied NE-SW or SW-NE clusters (for simphcity we draw them as 
straight hues, recall however that they will in general bend since they have 
two bending directions). Each of these clusters is blocked since it ends in a T- 
junction with a cluster along the transverse diagonal. The dotted rectangle 
adjacent to cluster AB (EF) are the regions in which this cluster can be 
displaced and yet a frozen backbone is preserved. Indeed the T-junctions in 
C and D (G and H) will be displaced but none of them will be disrupted. 
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can be expanded until emptying the whole lattice. From the results in Sec- 
tion 16.21 it is easy to see that the probability to expand an empty nucleus 
to infinity is dominated by the probability of expanding it up to n = ^y. 
Indeed, above this size the probability of an event which prevents expansion 
is exponentially suppressed. Therefore, considering the 0{L/^\\Y possible 
position for a regions that it is guaranteed to be emptyable up to size ^y, we 
get that the probability that is internally spanned is roughly bounded 
as R{L^ p) > L^S for < 6, where 6 is the probability that a small empty 
nucleus can be expanded until size ^y. In the emptying procedure described 
in Section we subsequently increase from n to n + 2 the size of the empty 
nucleus. If we require instead that on a region of size around each 

corner there is not a spanning DP cluster, we can expand directly to ^y. The 
probability of the latter event gives 6 > exp(— o^||~^) (see [15]). Therefore, 
we get the desired upper bound on logS that has the same scaling as the 
lower bound. 

6.5 Numerical Results 

In the following we present numerical results on the percolation of blocked 
structures that support our theoretical findings. Starting from an initial 
configuration we sequentially pick away all particles that are not blocked, 
so that the final configuration contains the backbones of forever blocked 
particles. In Fig. ^Jwe show the probability that a site belongs to the infinite 
spanning cluster for lattice sizes N = 100^,200^,400^,800^,1600^ averaged 
over 40000 samples. Though the transition takes place at the critical density 
of DP Pfip ~ 0.701, as we have proved, finite size effects appear already at 
0.53 - 0.55. 

Indeed, the probability that there exists a frozen cluster is substantial 
for p fifteen percent below pc even in our largest systems, (L = 1600): it is 
thus hard to study the asymptotic critical behavior (see [58] for an analogous 
problem in the context of bootstrap percolation). But in a slightly different 
model one can get closer to the transition [33]: these data are consistent with 
the predicted InS ~ {pc — p)~^ with p = 0.64, but the small range of InL 
available makes the uncertainties in p large. 

Another theoretical finding confirmed by numerics is the first order char- 
acter of the transition. Indeed the histograms of the density of blocked 
structures, ph, clearly show a well defined two peak shape as in usual first 
order transitions (see inset of Fig. IT^ . The peak at a non zero density is 
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Figure 13: Probability that a site belongs to the infinite blocked span- 
ning cluster as a function of the initial density for system sizes = 
100^2002,400^800^16002 averaged over 40000 initial configurations (from 
left to right respectively). Inset: histograms of the number of initial configu- 
rations (y axis) leading to a fraction of sites (x-axis) for = 800^, 40000 
initial configurations and initial density p = 0.5775, 0.5825, 0.5875, 0.5925 
(when increasing the initial density the right peak moves to the right). 
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Figure 14: Density density correlation as a function of time (Mon- 
teCarlo steps) for a 100 x 100 lattice. From bottom to top, p = 
0.49, 0.50, 0.51, 0.52, 0.53, 0.54, 0.55, 0.56, 0.57. 

clearly distinct from the one at zero density and moves to the right by in- 
creasing the density of particles in the initial configuration, thus showing that 
asymptotically the transition is first order. Furthermore there are no finite 
size effects on the peak position but only on its weight. Thus, the continuous 
curves Pb{p) in Fig. El are the average of two curves: the trivial one, Pb = 0, 
corresponding to the left peak and the dotted one corresponding to the right 
peak. By changing the system size the weights of the two peaks change and 
so the continuous averaged Pb{p) shifts to the right. 

Finally, let's focus on the predictions for the dynamics. First, as already 
discussed, an ergodicity breaking transition occurs at p = pdp. Furthermore, 
the first-order character of the percolation of blocked structures implies a 
discontinuous jump of qEA{p) = limt-+oo < Vx(t)rix{0) >c- This is the plateau 
of the correlation function, see Fig^^ and is the analog of the Edwards- 
Anderson parameter in spin-glasses. In fact for p < p^, ergodicity implies 
that qEA{p) = since the system always relaxes to equilibrium. Instead 
for p > pc, by separating the contribution from sites which are occupied and 
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blocked forever and the remaining sites, one can show that qEA{p) > c{p)(l{p) 
[15,37], where c(p) is a strictly positive for all p < 1 and q{p), as defined 
previously, is the probability the the origin is frozen. Therefore, qEA{p) > 0, 
for p> pdp. 

Furthermore the relaxation timescale, r, diverges also very fast at the 
transition as discussed previously. As a check, we performed numerical sim- 
ulations of standard Montecarlo dynamics. The results for < rix[t)r]x{'^) >c 
are plotted in Fig. EJfor system size = 100^. The curves clearly shows 
a developing (discontinuous) plateau which becomes infinite after a certain 
(size dependent) density. However, the initial curves until p = 0.53 have no 
finite size effects. This shows that the developing of a discontinuous plateau 
and the increasing of the timescale start far from the critical density, in 
agreement with an essential singularity for the laws for S, r (for a power 
law divergence r should diverge much closer to pc). Note that, because of 
the very severe finite size effects, one needs much larger system sizes (and 
much larger timescales) in order to measure successfully the exponent of the 
essential singularity. We leave that for future study. 

7 Conclusion 

As discussed in this work cooperative KCMs display a remarkable physical 
behavior. In particular. Knights model undergoes a purely dynamical phase 
transition due to a jamming percolation : a giant blocked cluster appears 
at Pc < 1- This new type of transition has features that are very different 
from usual (first and second order) phase transitions and that are similar 
to the ones indeed expected for glass-jamming transitions. In particular 
the fraction of jammed particles is discontinuous (as for first order phase 
transition) although time and lengthscales diverge (as for second order phase 
transitions). Furthermore the relaxation time diverges with a Vogel-Fulcher 
like form, i.e. much faster than a power law. 

The extension and the universality of our results are fundamental open 
questions. In three dimensions, two natural generalizations of our jamming 
percolation exist: one composed of DP clusters — which should slow down 
as a double exponential of (p — Pc)~^ — and the other of directed sheet-like 
structures which will have exponential slowing down like we have found in 
2D. The key ingredients are kinetic constraints that enable huge jammed 
clusters to form out of small objects without these becoming much more 
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common or much larger. A very important issue is whether the mechanism 
that we devise to create a jamming percolation, which is based on interacting 
DP clusters, is the only possible one or there are others that may lead to a 
very similar dynamical transition, see [54]. 

For the future, the connection between our results and the jamming tran- 
sition found for continuum particle systems [2] needs exploring. One should 
analyze the effects of constraint-violating processes occurring with a very 
low rate: these are certainly present for molecular liquids undergoing a glass 
transition. Similarly, it would be very important to generalize our results to 
systems of particles in the continuum. Finally, it would be very interesting 
to compare the geometrical and statistical properties of jamming percolation 
clusters to the ones that can be measured in experiments in colloidal and 
granular systems. 
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